Structural distortion and the spin liquid state in Tb 2 Ti 2 07 
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It is shown that a k — 0, A2 U distortion of the terbium tetrahedral network in Tb2Ti207 causes 
the apparent isolation of single tetrahedra as seen in neutron scattering studies. Single tetrahedron 
collective spin states, rather than individual spins, account for the main features of the spin liquid 
state, namely, fluctuating local moments and the absence of long range order. Singlet and doublet 
collective spin ground states are considered. An effective interaction between tetrahedra on the fee 
lattice is derived and found to be weak and anisotropic. 
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Tb2Ti207 is an example of three dimensional geomet- 
ric frustration because the magnetic Tb 3+ ions are ar- 
ranged on a corner-sharing tetrahedral network. It has 
a spin liquid phase, characterised by the absence of long 
range correlations and rapidly fluctuating local magnetic 
moments, which persists down to at least 50 mKji^ 
This article is devoted to another mysterious feature of 
Tb2Ti207, namely the apparent isolation of single tetra- 
hedra from the tetrahedral network. Such isolation is 
suggested by neutron scattering experiments, which find 
that magnetic correlations beyond the size of a single 
tetrahedron are absent ?i and exact calculations of spin 
eigenstates on single tetrahedra reproduce well diffuse 
neutron scattering patterns.—*^ 

A reasonable description of the exchange interaction 
across the four Tb sites on a single tetrahedron in 
Tb2Ti2C>7 has been obtained only recently.-^ The crystal 
electric field (CEF) ground state of the Tb ions is a dou- 
blet, giving rise to a classical picture of Ising-like spins 
constrained to point into or out of the four vertices of 
the tetrahedron. However, mixing with higher CEF lev- 
els tends to restore transverse spin components, negating 
the classical picture. The resulting ground state resem- 
bles neither the ferromagnetic "spin ice" states, which 
have two spins pointing inside the tetrahedron and two 
pointing out, nor the anti- ferromagnetic "all-in or all- 
out" states. Instead it is a linear combination of various 
states and with no classical analogue. 

The issue of how to deal with the entire tetrahedral 
network, and why the single tetrahedron picture is valid, 
remains an outstanding problem. Because of the connec- 
tivity of the tetrahedral network, eigenstates of a single 
tetrahedron are not in general eigenstates of the entire 
tetrahedral network. This article offers a solution to this 
paradox, which involves symmetry breaking in the form 
of a lattice distortion. 

Lattice distortions are often invoked in theory as a 
way to relieve geometrical frustration and arrive at an 
ordered state,— and usually these involve a change of the 
crystal system, for example, from cubic to tetragonal. In 
fact, there is evidence that such a transition occurs in 
Tb2Ti2C>7 below any accessible temperature.— In this ar- 
ticle, we propose another kind of lattice distortion, one 



which lowers the point group symmetry but leaves the 
crystal system unchanged. This distortion does not re- 
move the frustration completely, but it does serve to iso- 
late single tetrahedra. Collective spin states on single 
tetrahedra then replace individual spins as the funda- 
mental basis for excitations and longer range effective 
interactions. The main features of the spin liquid phase 
follow from this scenario. 

Pyrochlore crystals such as Tb2Ti207 belong to the 
cubic space group Fd3m (Oh 7 , No. 227), in which both 
Tb and Ti ions form separate, corner-sharing tetrahedral 
networks. Here we are only concerned with the mag- 
netic Tb 3+ ions. The tetrahedra appear in two different 
orientations (A and B), which alternate in the tetrahe- 
dral network. A tetrahedron of either type comprises the 
primitive unit cell. Thus the set of all A tetrahedra form 
a face-centred cubic (fee) lattice, as do the set of all B 
tetrahedra. The set of A tetrahedra are related to the 
set of B tetrahedra by the 7r and 7r/2 screw rotation and 
inversion elements of FdZm. 

We begin by considering the nearest neighbour 
isotropic (Heisenberg) exchange interaction. The Hamil- 
tonian for the entire tetrahedral network can be split into 
two terms 



He x — J Ha + J Hb 
where J is the exchange coupling constant, 



H A =Y J X, 



(Ak) 



(1) 



(2) 



and X(Ak) is the exchange interaction over the k-th A 
tetrahedron, and Hb is the exchange interaction over the 
set of B tetrahedra. In either case, the sum over A; is a 
sum over fee lattice sites. The exchange interaction is 



(3) 



<i]> 



where the sum runs over four sites on a given tetrahe- 
dron. The eigenstates of Ha and Hb are known, and are 
simply the direct products across tetrahedra of the eigen- 
states for a single tetrahedron, which are given in Ref. 4 
and described below. However, in general, [Ha, Hb] ^ 0, 
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therefore Ha and Hb do not have common eigenstates 
with each other or with H ex (an exception to this is dis- 
cussed below) . This is in apparent contradiction with the 
results of experiments, which agree very well with the sin- 
gle tetrahedron picturei^ This suggests that H ex dJ) 
should be replaced by an effective Hamiltonian^ 

H eff =J A H A + J B H B (4) 

where one of the couplings is much larger than the other. 
The difference between the exchange couplings J a and 
Jb could originate from a structural distortion which 
makes the A tetrahedra smaller and the B tetrahedra 
larger, thus J a > Jb, or vice versa. Such a distor- 
tion is a k — 0, A 2u mode of the Tb ions which reduces 
the space group symmetry from Fd3m to -F43m.— This 
new space group is compatible with an otherwise unex- 
plained (0, 0, 2) Bragg peak observed in neutron scatter- 
ing measurements^ If present, domains could produce 
isolated spins along their boundaries, possibly contribut- 
ing to observed glassy behaviour ^ 11 ' 12 

Fits to the dominant (3,1,1), (2,2,2), (1,1,1) and 
(2, 2, 0) Bragg peaks are in good agreement with the py- 
rochlore space group Fd3m^ Therefore it is reasonable 
to assume that the much smaller (0,0,2) peak is due to 
a small distortion of Fd3m. The space groups which are 
compatible with all five of the peaks, are cubic, and which 
are also subgroups of Fd3m are F43m (No. 216, T 2 ,), 
P43m (No. 215, T d ), F23 (No. 196, T 2 ), P23 (No. 195, 
T 1 ) and P2i3 (No. 198, T 4 ). Among these, F43m has the 
highest symmetry, and all are subgroups of FAim. Thus 
it is strongly indicated that the space group of Tb2Ti2C>7 
is in fact F43m. The k = 0, A 2u displacement mode 
which gives rise to F43m can occur for the Tb, Ti or 
O ions,— although physical considerations point in favour 
of Tb displacements, as discussed above, or O displace- 
ments due to the role that the O ions play in mediating 
the exchange interaction. Whether or not the (0, 0, 2) 
peak is present at all temperatures, or appears due to 
some kind of exotic "spin Jahn- Teller" effect, remains to 
be investigated^ 

The ratio Ja/Jb varies with the size of the lattice 
distortion. The localised nature of the Tb 4/ electrons 
ensures that overlap integrals contributing to the ex- 
change constants will be very sensitive to changes in 
relative distances, and this is evidenced by very large 
magnetostriction.— Therefore it is possible that even a 
very small lattice distortion could produce a large differ- 
ence between J a and Jb- 

Assuming that \Ja\ > \Jb\, we will find the ground 
state of Ha and then consider Hb as a perturbation. 
In order to do this, we must first describe the collective 
spin states of a single tetrahedron. A tetrahedron has 
four magnetic ions at its vertices, each of which has a 
local site symmetry of D^d- Following the conventions 
established in Ref. 4, we label the ions #1, 2, 3 and 
4, where their C3 axes point in the directions [1,1,1], 
[-1, -1, 1], [-1, 1, -1] and [1, -1, -1] respectively. Tb 3 + 



ions have a total angular momentum of J — 6, but the 13- 
fold degeneracy is split by the crystal electric field into 5 
singlets and 4 doublets; one of the doublets is the ground 
state. We take the ground state asi^ |±) = ±0.13| ± 5) =F 
0.13] =F 1) — 0.95 1 T4), where the quantisation axis points 
in the direction of the C3 axis for each Tb ion. There 
are then 16 collective spin states on a tetrahedron, which 
can be written as | ± ±± ±) = |±)i ® |±) 2 ® |±) 3 ® |±) 4 , 
where the subscripts indicate the site on the tetrahedron. 
Symmetry considerations predict that interactions will 
split the 16 states into a singlet, three doublets and three 
triplets, 4 which are complicated linear combinations of 
the basis states | ± ± ± ±). 

The ground state spin configuration of Tb ions on a sin- 
gle tetrahedron for isotropic anti-ferromagnetic exchange 
is a doublet, 4 which we write as \E±). Otherwise the 
ground state depends on the details of the anisotropy. 
An anisotropic interaction (equivalent to including near- 
est neighbour dipole-dipole interactions) was considered 
in Ref. 5! and the singlet was found to be the ground state. 
A triplet ground state could also yield patterns similar to 
what are found in experiment A In the following, we will 
consider the doublet and the singlet as possible ground 
states on a single tetrahedron. 

The ground state of Ha is constructed by taking the 
direct product of the single tetrahedron ground states 
across the A tetrahedra. If the tetrahedron ground state 
is the singlet, 

l^> = (! + + --> + 1 -- + +) + ! + - + -) 

+ |- + -+> + l + - -+> + |- + + -»A/6 (5) 

then the ground state of Ha can be written 14^1^!)^, 
where k indexes tetrahedra, that is, fee lattice sites. This 
state clearly has long range correlations but they will be 
undetectable in static neutron scattering measurements, 
in agreement with experiment.— This is because neu- 
tron scattering detects magnetic correlations, which are 
proportional to the matrix elements (^l Ji|$)( ( E ) | Jj|\E r ), 
where i and j are Tb sites and ^ and <E> are eigenstates 
of Ha, which vanish unless i and j belong to the same 
A tetrahedron. On the other hand, higher order corre- 
lations, beginning with quadrupolar, will exist between 
tetrahedra. Thus the state IIfc|Ai)fe is not a true spin 
liquid. The excitation spectrum of Ha will be gapped, 
with an energy of the order of 0.1 ~ 1 corresponding 
to the separation between |^4i) and the first excited state 
on the tetrahedron. 

If the tetrahedron ground state is a doublet \E±) then 
the ground states of Ha can be written as Ylk \Eo)k- 
The ground state is highly degenerate. One way of view- 
ing the degeneracy is to note that the state at each site 
can be any complex linear combination of the doublet 
states \E±). This SU(2) freedom on each fee site leads to 
the absence of correlations of any kind beyond the size 
of a single tetrahedron. However, interactions coming 
from the perturbation Hb will limit this freedom. Ulti- 
mately, we find a system of weakly interacting tetrahedra 
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arranged on an fee lattice (which is also frustrated). If 
the tetrahedron ground state is a triplet then there will 
be a SU(3) symmetry on each lattice site, which could 
also be limited by weak interactions. 

In order to calculate the effect of the perturbation Hb 
we first need to examine in detail the exchange inter- 
action over a single tetrahedron and its eigenstates. A 
useful representation of the exchange interaction over a 
single tetrahedron is given in Ref. 3 in which the angu- 
lar momentum operators are expressed in terms of local 
coordinate axes (indicated by subscripts) such that the 
local z axis for each Tb ion points in the direction of 
its C3 axis. Local x and y axes have also been implicitly 
selected. The result is written in Table I. Only local coor- 
dinate operators are used in the following. The symmetry 
of the crystal is octahedral, which permits three separate 
invariant terms in the exchange interaction, J\X\, J 2 X 2 
and J73X3, where c7i,2,3 are different exchange couplings 
for each term, and X±^,s are the sums of all the terms 
in the second, third and fourth columns respectively of 
Table I. There is actually a fourth invariant allowed un- 
der Oh symmetry (A3 is split), which we do not consider 
here, although the following discussion can be easily gen- 
eralised. The isotropic exchange interaction is the sum 
of all three terms with J\ — J 2 — J?,- In the follow- 
ing, we will assume that the exchange interaction is anti- 
ferromagnetic. 

The matrix elements for J z ,± arc 

(±\J Z \±) = ±j (6) 
(±\J±\T) = t (7) 

The parameters j and t serve as a very useful characteri- 
sation of single-ion doublet spin states. For J = 1/2 ions, 
j = 1/2 and t = 1, while for other 1/2-integral spins j 
and t may be quite different. The parameter t equals zero 
for integral spins (as for the Tb 3+ ion) but may acquire a 
significant non-zero value due to mixing with higher crys- 
tal electric field levels^ A comparison between theoreti- 
cal and experimental diffuse neutron scattering patterns 
suggests that t is larger than j in Tb2Ti2C>7. 4 ' 5 Thus, 
in our initial approximation, we will assume that j = 0; 
then only term in the exchange interaction with non-zero 
matrix elements is A3. 

The ground state of A3 is the doublet^^ 

\E ± ) = V375|4 1} > - >/W£ 3) > (8) 

where 

l4 1} > = ! + + + +>, \e w ) = \ ) (9) 

|4 3) > = (I + + - -)+e\ + - + -)+s 2 \+- -+) 
+ | -- + +)+ e\ - + -+)+ e 2 \ - + + -) )/V6 (10) 

|£(3) } = [|s (3) )r (n) 

where e = exp =|p First order corrections due to Hb are 
found by calculating the matrix elements of Hb between 



the degenerate ground states of Ha- We will express 
the final result as an effective interaction between the 
ground state doublets on neighbouring tetrahedra. Re- 
call that Hb is the exchange interaction summed over all 
B tetrahedra. It consists of terms of the form J; • Jj, 
where i and j are nearest neighbour but are found on 
different, neighbouring A tetrahedra. Nearest neighbour 
A tetrahedra share exactly one pair of nearest neighbour 
magnetic ions. Considering only the A3 part of the ex- 
change interaction, we find that all of the ground state 
matrix elements vanish^ because each term in Hb raises 
or lowers at most one site in any A tetrahedron, and the 
result is orthogonal to all of the ground states. There- 
fore, to first order in JbI Ja and zeroth order in j/t, the 
A tetrahedra are non-interacting. 

Now we consider finite values of j/t by including Xi 
and A2 in the exchange interaction. This will add mix- 
tures of 

l4 2) > - (1 + — ) + !- + --) 

+1 -- + -> + I +»/2 (12) 

|4 2) ) = (| - + + +) + 1 + - + ++) 

+| + + -+) + I + + + -»/2 (13) 

and 

\Ep) = VV5\E { ± ] ) + VySlE^} (14) 

to \E ± ). The general form is \E±) = a\E { ± ] ) + /3\E^) + 

5\E± ). To find the effective interaction between tetra- 
hedra, consider two neighbouring A tetrahedra, with Tb 
sites numbering 1, 2, 3, and 4 on the first tetrahedron 
and 5, 6, 7 and 8 on the second. Projecting all possi- 
ble ground states of Ha onto the subspace of these two 
tetrahedra yields four states of the form \E±; E±). Sup- 
pose that sites 1 and 6 are nearest neighbours and that 
the C3 axis of site 6 points in the same direction as that 
of site 2. Then J\ ■ Jg is a term in Hb and takes the 
same form as J± ■ J2- The effective interaction between 
tetrahedra is found by calculating the 16 matrix elements 
(E±; E±\Ji ■ Je\E±; E±) . The non-zero matrix elements 
are 

{E+;E±\X 1 \E+;E ± ) = (£_; E T \X X \E-; E T ) = T 7 2 /3 
(E+;E ± \X 2 \E +] E T ) = (E±; E+\X 2 \E T ; E+) = - 7 rV2/3 
{E_-E±\X 2 \E^-E T ) = (E ± ;E_\X 2 \E T ;E_) = 7 rV2/3 
(E + ;E + \X 3 \E_;E_) = (£_; E_\X 3 \E + ; E+) = r 2 /3 
(E + : E_ I A 3 1 E_;E + ) = (£_; E + \X 3 \E + ; E_) = -t 2 /6 

where 7 = j(a 2 — f3 2 /2) and r ~ ta(3. From these ma- 
trix elements we can infer that the effective interaction 
between tetrahedra is anisotropic in general and can be 
divided into three separate terms 

Htetra = 3 \X\ + J ' 2 X 2 + J' '3X3, (15) 

where A, now operate on tetrahedral ground state dou- 
blets, as follows. There are four invariant operators on 



4 



term 


Xi 


x 2 
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J 2 
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Ji 
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- +E 2 Jl-)Jiz] 


^(e 2 Ji+J 3 + +eJi_J 3 - 
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-)• 
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TABLE I: The exchange interaction over a single tetrahedron expressed in terms of local coordinates for each Tb ion. The 
first column lists the terms in the exchange interaction over a single tetrahedron. In each row, the first entry is the sum of the 
middle three. The isotropic exchange is the sum of all the terms in the middle three columns, and may be divided into three 
anisotropic terms Xi, X 2 and A3 which are sum of all the terms in each of the middle three columns. The last column lists 
the bond direction for each term, e = exp 2hl 



the fee lattice that are directly related to the four invari- 
ants on the single tetrahedron (note that, as explained 
above, only three have been considered here). The bond 
direction between a tetrahedron and each of its twelve 
nearest neighbours can be defined by the bond direction 
between the nearest neighbour ion sites on neighbour- 
ing tetrahedra. Six different bond directions are present. 
The bond direction determines the correspondence be- 
tween rows in Table I and terms in the effective interac- 
tion between tetrahedra. Assuming that matrix elements 
of the tetrahedron operators are {E±\J Z \E±) = ±1 and 
(E±\J±\E T ) = 1, we find J\ = J l7 2 , J' 2 = J 2l T and 
J' 3 = J3T 2 . Higher order corrections arising from mix- 



ing via Hb between the single tetrahedron ground states 
\E±) and excited states will further renormalise J' \, J'2 
and J' '3. The renormalisation of the coupling constants 
in Eq. [T5] implies that even if the underlying exchange in- 
teraction is isotropic, the effective exchange interaction 
between tetrahedra is anisotropic. In a similar fashion, 
the effective interaction between tetrahedra with triply 
degenerate ground states can also be found. 

Anisotropy in H te tra may assist long range ordering 
of the tetrahedra. However, at temperatures so far ob- 
tained, there is no evidence from neutron scattering for 
any kind of inter-tetrahedra correlations. The inelastic 
neutron scattering function is proportional tc 19 



J(q)c< 



,-E m /k 



q-Ob-rv) 



(16) 



i,j a,b n 



where a and b are the four magnetic ion sites at the corners of a tetrahedron and n and m are eigenstates of the 
system. The angular momentum operators J 1 refer to the global coordinate system (indicated by superscripts). This 
formula yields patterns which are linear combinations of the following functions 



f(h,k,l) = 
f(h,k,l) = 
f(h,k,l) = 
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h 2 - 


h k 2 + 1 2 
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h 2 - 


Yh 2 +P 




1 



h 2 + k 2 + 1 2 



. Trh . nk . 7rfc . ttI irh . irl 

hk sm — sm 1- kl sm — sm V hi sm — sm — 

2 2 2 2 2 2 

9 7rfc irl 2 nil 7rZ 9 irh nk 

h cos — cos Vk cos — cos Yr cos — cos — 

2 2 2 2 2 2 



(17) 
(18) 



nh 



irk 



nk 



nil 



nl 



h cos — cos h cos — + k cos COS h cos — + I cos — cos h cos ( 



nl 



nh 



Inter-tetrahedron correlations would be manifested in the 
appearance of smaller period (in k space) contributions 
to the scattering patterns. If present, these could help 
to constrain the effective interaction between tetrahedra. 
However, their absence indicates that inter-tetrahedra in- 
teractions are indeed weak. 

Finally, let us contrast our results to situations in 
which the commutator [Ha,Hb] does vanish, allowing 
common eigenstates of Ha, Hb and H ex . Non-zero terms 



in the commutator [Ha,Hb] arise from terms propor- 
tional to X 2 and X 3 in Ha and Hb- Thus, if the exchange 
interaction is highly anisotropic and only X\ appears or 
when the parameter t vanishes, then [Ha,Hb] — 0, and 
the eigenstates of H, Ha or Hb are the sixteen basis 
states I ± ± ± ±). If the sign of J\ is positive then the 
ground states of a single tetrahedron are the antiferro- 

magnetic states | + + + +) and | ); otherwise 

the ground state has a six fold degeneracy, | + H ), 
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I — + +>» I + - + ->, I - + -+>, I + and | - + + -). 

These situations can each be defined by rules: "all in or 
all out" in the former and "two in/two out" (spin ice 
rule) in the latter. Ground states of Ha are again found 
by taking the direct product of ground states over the A 
tetrahedra. A ground state of Ha will be an eigenstate of 
Hb but in general it will not be a ground state unless the 
appropriate ground state rule is satisfied on all of the B 
tetrahedra too. The distinguishing feature between the 
special case when [Ha,Hb] = and the general case is 
the form that the eigenstates take. When the commuta- 
tor is non-zero the eigenstates of the single tetrahedron 
are necessarily entangled, that is, they must be linear 
combinations of the basis states | ± ± ± ±) , which leads 
to fluctuating local moments found in Tb2Ti207. 

To summarise, we have shown that a k = 0, A2 U lattice 
distortion can account for the observed single tetrahedron 
behaviour. This distortion results in a space group of 



F43m, which is compatible with with a (0, 0, 2) Bragg 
peak observed by neutron scattering^ and the loss of 
inversion centre recently observed in Raman scattering.— 
The eigenstates of the system are then direct products 
over the fee lattice of the single tetrahedron eigenstates. 
The effective interaction between tetrahedra is weak and 
anisotropic. 
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